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Abstract 

The new measurements of B — > J/tpK(K*) decays by CDF and CLEO indicate that 
the production ratio R and the fraction of longitudinal polarization Tl/T are smaller 
than the previous results. In conjunction with the new result of parity-odd transverse 
polarization in B —* J/^K*, we found a minimal modification to the factorization 
hypothesis: While the data of B — > J '/iftK* can be accommodated in the factorization 
approach with nonfactorizable terms x A = Xa 2 = X v — X of order 15%, the result of 
R measurement requires that the nonfactorizable effect x Fl on B — > J/^K be slightly 
larger than \. Therefore, the effective parameter a| is not universal even for B — ► 
J/ifjK(K*) decays. We have generalized the considerations to B — > ip(2S)K(K*) and 
B s — > J/ifxj) and found that the predictions are in agreement with currently available 
data. 
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The two-body nonleptonic weak decays of D and B mesons are conventionally described 
by the factorization approach. It has been shown [|], |2| that this approach fails to account 
for the observed fraction of longitudinal polarization Tl/T in B — > J/ipK* decays in all 
commonly used models of form factors and the data of the production ratio R = T(B — > 
J/ijjK*)/T(B — > J/ipK) in many known models. The issue of how to test or modify the 
factorization hypothesis to describe the data of B — > J/ipK(K*) has been the subject of many 
subsequent studies. Recently, CDF || and CLEO [|J have presented new measurements of 
B — > J/ipK(K*) decays and for the first time CLEO has measured the parity-odd transverse 
polarization in B — > J/ipK* . It turns out that the new results of Tl/T and R are smaller 
than the previous values. The purpose of this Letter is to study the theoretical implications 
of the new data. 

Under the factorization approximation, the hadronic matrix element is factorized into 
the product of two matrix elements of single currents, governed by decay constants and form 
factors. For the QCD-corrected weak Hamiltonian 

H cS oc ciOi + c 2 2 = ci(gig 2 )(g 3 g 4 ) + c 2 (gig 4 )(g3<?2), (1) 

where (qiq 2 ) = <?i7^(l — 75)92, the external W-emission and internal W^-emission amplitudes 
(or the so-called class-I and class-II decay amplitudes) characterized by the parameters a>i, a 2 
are related to the Wilson coefficient functions c\ and c 2 via 

a 1 = c 1 + -, a 2 = c 2 + -. (2) 

Factorization requires that a\ and a 2 be universal; i.e., they are channel-independent in D 
or B decays. However, we have learned from charm decays a long time ago that the naive 
factorization approach never works for the decay rate of color-suppressed (i.e., class-II) decay 
modes, though it might work for class-I decays f| . The most noticeable example is the ratio 
T(D° — > K°tt°)/T(D — > K~tt + ), which is predicted to be ~ 0.02 whereas experimentally 
it is measured to be 0.55 ± 0.06 ||. This implies that the inclusion of nonfactorizable 
contributions is inevitable and necessary. Since the amplitudes of D, B — > PP, VP decays 
(P: pseudoscalar meson, V: vecetor meson) are governed by a single form factor, the effects 
of nonfactorization amount to a redefintion of the effective parameters a\ and a 2 0: 

af = Cl + c 2 ^ + Xi) , af = c 2 + Cl (— + x 2 ) , (3) 

or 

af = a x (l + ^ X i) , af = a 2 (l + ^2) . (4) 

For the expressions of the nonfactorizable terms denoted by the parameters xi an d X2, 
consider the decay B — > Dn as an example; we obtain || [| 

v ,s n x v (g^rjQ^B) fll (DttIOxIB)^ F (8)w/ (m^) 0l if^K) 



^o SD Ptt|0 1 | j B) / c 2 ( J D7r|0 1 | J B} / F ^«) c 2 F BD (ml] 
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v ,5 n] = v (^vr|d 2 |g) a 2 (D7r|0 2 |B) B/ _ ffj^mgQ a 2 ^V|) 

Xzl J (£>7r|0 2 | J B} / + Cl (DTrldlS), F^(ml,) + ci Ff*{m%) ' W 

where the subscripts / and nf denote factorizable and nonfactorizable contributions, respec- 
tively, and use has been made of the Fierz identify 

Oi,2 = j0 2l i + 2)1 (6) 

with Oi = |(<?iA a g 2 )(g 3 A a g 4 ) being a product of the color-octet currents (q\ a q') = qj^l — 
75)A a g'. Since the factorizable B — »■ Z^7r amplitudes are 



(£>7r|0 1 | J B) / = ^(m^-m^Fo^K), 

(IhrlOalB), = f D (m 2 B -ml)Ft(m 2 D ), (7) 
we have followed [[J to define the form factors F ^ and Fq 8 ^^: 

(2>7r|Oi|S>^ = /.K-m^F^K), 
(2>7r|Oi|S> = /.(m|-m^)F (8)n/ (m*)- (8) 
Likewise for i? — > Dp and _D*7r decays, we have 

-»• Dp) = x ifBD , x 2 (B^Dp) = X Bp , 
Xl (S - D*ir) = X , x 2 (B -> -D*tt) = X 2 fS „ . (9) 

Since ci/a 2 >> 1, it is clear from Eqs.(4) and (5) that nonfactorizable effects are more 
dramatic in color-suppressed decay modes. As a consequence, a determination of X2 is more 
reliable than xi- 

The study of nonfactorizable effects in M — > V V decay is more complicated as its general 
amplitude consists of three independent Lorentz scalars: 

A[M{p) - Pi(£i,pi)y 2 (£ 2) p 2 )] oc£;(A 1 )£:(A 2 )(i 1 ^ i/ + i 2 pV + ^e Aim Va^), (10) 

where Ax> A 2 , V are related to the form factors Ai, A 2 and V, corresponding to S-, D- and 
P-waves, respectively. For B — > D*p we can define various nonfactorizable terms: 



X 1 ,. 2 (5^ J D»=X 1AflD ,, X 2 , A2 (B^D*p) = X 2ABp , (11) 

Since a priori there is no reason to expect that nonfactorizable terms weight in the same way 
to S-, P- and D-waves, namely x iAl = X iA2 = X iV i ^ * s * n general not possible to define 
an effective a\ or a 2 for M — > VV decays once nonfactorizable effects are taken into account 



T0|1 • However, if nonfactorizable contributions are channel-independent, i.e., 

Xi,F Xi l F 1 Xi,A Xi l A 1 Xi,A 2 "X.i.y X? 

(12) 
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the effective parameters af 1 and a| ff given by Eq.(3) or (4) are applicable not only to 
M — > PP, VP but also to M -> VV and they become universal. In this case we have 
a new factorization scheme with x ^ 0. Note that the predictions for the ratios of physical 
quantities (e.g., Y^/Y and |P| 2 ) in the same class of decay modes in the new factorization 
method are the same as that in naive factorization (i.e., % = 0). Empirically, it has been 
found that that the discrepancy between theory and experiment for charm decays is greatly 
improved if Fierz transformed terms in (3) are dropped |5|. It has been argued that 
this empirical observation is justified in the so-called large- N c approach in which a rule of 
discarding sub leading 1/N C terms can be formulated jL^]. This amounts to having universal 
nonfactorizable terms xi = X2 = —^/N c in Eq.(3). In reality, there is no reason to expect a 
universal Xi or X2- Since xi comes mainly from color-octet currents, it is natural to expect 
that 0, 



\ X 2{D - VP) I > \ X 2{D - PP)| > \ X *(B - VP) I, 



(13) 



as nonfactorizable soft-gluon effects become more important when final particles move slower, 
allowing more time for significant final-state interactions. Indeed, the above theoretical 
expectation is confirmed by the phenomenological analyses of D and B decay data 0, [T^ . 
For charm decay, Xi is found to be in the range —0.60 < x%(D) < —1/3 [fl3[1 . 

We now come back to B — > J/ipK(K*) decays. The general expressions for the fraction 
of longitudinal polarization Yl/Y and production ratio R read 



Yl 
Y 

R 



Y(B -> J/iPK*) l 
Y(B — > J/i)K*) 
Y(B -> J/^K*) 
Y(B -> J/ipK) = 



a \1 
1.08 



Cl 

a 2 



-X; 



( 1+ 



bx[l + —x Ar 
a 2 ' 



2 ! 



(14) 



where 



x 



" ( 1 + — Xa, 

a 2 1 



A 



BK* 



,171 



Jli>> 



bx(l + —Xj 

\ a 2 

V BK* {m 2 



Cl 

1 + a 2 ^ 



= Ft 



1 2 2 
+ cy 



BK( m%^) 



1 + ^X V 
a 2 



(15) 



and Xa 1 denotes X 2A bk* e ^ c -' f° r example [see Eq.(5)], 

4 (8)n// 2 \ 4 (0)n/. 

A i l m jA/J a 2 A l ( 



X., 



X., 



l(8)n vw 



+ 



m 



j/i>) 



m 



(16) 



The analytic expressions of a, b and c are given in [III El. Numerically, 

a = 3.164, 6= 1.304, c = 0.435 . 



(17) 
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The parity-odd (P-wave) transverse polarization measured in the transversity basis [14 
which is more suitable for parity analysis, has the form 

|P| 2 





A ± 


2 




2 + 




2 + 


A± 


2 



where A^s are defined by the orientation of the J/ip polarization vector ej/^ in the transver- 
sity basis ||]: A for ej/y, parallel to x, A\\ for ej/y, parallel to y, and Aj_ for ej/^, parallel to 
z. Since A^'s are related to the amplitudes H\ in the helicity basis via 

H ± = ^=(A ll ±A ± ), H = -A , (19) 



P\ 2 =2c 2 y 2 (l + ^ Xv ) It (20) 



we find 

X.v 

a 2 

A measurement of |P| 2 will thus provide information on the vector form factor. 

We would first like to see if the data of Tl/T, R and |P| 2 can be explained in the 
factorization approach. As we have stressed in passing, naive factorization does not work 
for the decay rate of color-suppressed decay modes; nonfactorizable contributions should 
always be included in order to describe the branching ratios of B — > J/ipK(K*). However, if 
Xf 1 = X Al = Xa 2 = Xvi we nave a new factorization scheme but the predictions of Tl/T, R 
and \P\ 2 in the naive factorization method remains intact since all nonfactorizable terms 
are canceled out in Eqs.(14) and (20). To proceed, we consider several phenomenological 
models of form factors: (1) the Bauer-Stech-Wirbel model (BSWI) |15[ in which B — ► 
K(K*) form factors are first evaluated at q 2 = and then extrapolated to finite q 2 using a 
monopole behavior for all form factors, (2) the modified BSW model (BSWII) |Tj|, which 
is the same as BSWI except for a dipole q 2 dependence for form factors Fi, A , A 2 and 
V, (3) the nonrelativistic quark model by Isgur, Scora, Grinstein and Wise (ISGW) |L7| 
with exponential q 2 dependence for all form factors, (4) the model of Casalbuoni et al. and 
Deandrea et al. (CDDFGN) |18| in which form factors are first evaluated at q 2 = using 
heavy meson effective chiral Lagrangians, which incorporate heavy mesons, light pseudoscalar 
mesons and light vector mesons, and then extrapolated with monopole behavior. Several 
authors have derived the B — ► K(K*) form factors from experimentally measured D — > 
K(K*) form factors at q 2 = using the Isgur- Wise scaling laws based on heavy quark 
symmetry [19], which are allowed to relate B and D form factors at q 2 near g^ ax . The 



B — > K(K*) form factors are calculated in Q [IW(i)] by assuming a monople dependence 
for all form factors, while they are computed in j2(J [IW(ii)] by advocating a monopole 
extrapolation for F±, Aq, A±, a dipole behavior for A 2 , V, and an approximately constant 
for Fq. An ansatz proposed in ||, which we call IW(iii), relies on "soft" Isgur- Wise scaling 
laws and different q 2 behavior of A\ from A 2 , V, Fx- 
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In all above form-factor models, the q 2 dependence of form factors is assumed to be 
governed by near pole (monopole or dipole) dominance. It is thus important to have a 
first-principles or model calculation of the form-factor q 2 dependence. In principle, QCD 
sum rules, lattice QCD simulations, and quark models allow one to compute form-factor 
q 2 behavior. However, the analyses of the QCD sum rule yield some contradicting results. 
For example, while A Bp is found to decrease from q 2 = to q 2 = 15 GeV 2 in [21], such a 



phenomenon is not seen in p2| , |23|j . Also the sum- rule results become less reliable at large 
q 2 due to a large cancellation between different terms. The present lattice QCD technique 
is not directly applicable to the B meson. Additional assumptions on extrapolation from 
charm to bottom scales and from g^ ax to other q 2 have to be made. As for the quark 
model, a consistent treatment of the relativistic effects of the quark motion and spin in a 
bound state is a main issue of the relativistic quark model. To our knowledge, the light- 
front quark model |24|] is the only relativistic quark model in which a consistent and fully 



relativistic treatment of quark spins and the center-of-mass motion can be carried out. A 
direct calculation of P — > P and P —>■ V form factors at time-like momentum transfer 
in the relativistic light-front (LF) quark model just became available recently p5] . The 
contributions from valence-quark configuration to B — > K(K*) form factors are p5| : Q 

Ai K * (m 2 ^) = 0.43, V BK * (m 2 ^) = 0.50. (21) 

We found that Ff K , Af K *, A BK * , V BK * exhibit a dipole behavior, while Af K * shows a 
monopole dependence in the close vicinity of q 2 = [P3J . Table I summerizes the predictions 
of Tl/T, R and \P\ 2 in above-mentioned various form-factor models within the factorization 
approach by assuming the absence of inelastic final-state interactions. 
The pre- 1996 values of R and Tl/T are given by 

R= 1.68 ±0.33, T L /r = 0.74 ±0.07, (22) 

where the former is the average value of ARGUS 2j| and CLEO [p7| , and the latter is the 



combined average of ARGUS: T L /T = 0.97±0.16±0.15 ||, CLEO: T L /T = 0.80±0.08±0.05 
27], and CDF: T L /T = 0.65 ± 0.10 ± 0.04 |2j§. It is obvious from Table I that all the 



existing models fail to produce a large longitudinal polarization fraction, whereas several 

^dn the light-front model calculations, the valence-quark contributions to form factors Aq, A%, V depend 
on the recoiling direction of the K{K*) relative to the B meson [p5| . Thus the inclusion of the non- valence 
configuration arising from quark-pair creation is in principle necessary in order to ensure that the physical 
form factors are independent of the recoiling direction. Since the non- valence contribution is most important 
only near zero recoil, we expect that Z-graph effects on B — > K(K*) form factors at q 2 — m 2 j/^ obtained in 
the "+" frame, where K* is moving in the +z direction in the rest frame of the B meson, are not important. 
As noted in p5| , the behavior of V BK in the "+" frame is peculiar in the sense that in general the form 
factor in the "+" frame is larger than that in the "— " frame, so we have taken V BK from the "— " frame 
and other form factors from the "+" frame. 
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Table I. Predictions for form-factor ratios x, y, z and for Tl/T, R and \P\ 2 in various form-factor 
models. 



BSWI p| 



BSWII JT6| 
ISGW [0 
CDDFGN 
IW(i) [0 



IW(ii) [0 

iw(iii) g 

LF pBl 



pre-1996 expt. 



CDF [[§, 
CLEO [f 



2/ 







1.01 
1.41 
2.00 
1.00 
0.98 
0.88 



1.19 1.23 

1.77 1.82 

2.55 2.30 

3.24 2.60 



2.56 
1.77 



1.08 2.16 
1.16 1.35 



1.74 
2.05 
1.86 
1.78 



r L /r 



0.57 
0.36 
0.07 
0.37 
0.45 
0.56 
0.45 
0.50 



R 



4.22 
1.63 
1.72 
1.50 
2.89 
1.84 
2.15 
1.84 



0.74 ±0.07 1.68 ±0.33 



\P\ 



0.09 
0.24 
0.52 
0.30 
0.31 
0.16 
0.26 
0.13 



0.65 ±0.11 1.32 ±0.28 

0.52 ±0.08 1.36 ±0.20 0.16 ±0.09 



models give satisfactory results for the production ratio. Since the prediction of Yl/Y in 
new factorization with \a = Xa 2 = X v * s the same as that in naive factorization, this means 
that if the longitudinal polarization fraction is as large as 0.74 ± 0.07, then nonfactorizable 
terms should contribute differently to S-, P- and D-wave amplitudes. Consequently, the 
effective parameter a| ff cannot be defined for B — ► J /ipK* decays if Yl/Y is large. In order to 
produce a large Yl/Y, various possibilities of nonfactorizable contributions to B — > J/ipK*, 
e.g., S'-wave dominance: Xa ^ ^> Xa 2 = X v = ®i have been explored in [10, [13], |29| . 



The new results from CDF [^J on the branching-ratio measurements of B — > J /%pK{K*) 
decays give rise to 

R= 1.32 ±0.23 ±0.16. (23) 

Recently, CLEO II |3J has completed the analysis of all data sample and come out with the 
results: 

R = 1.36 ±0.17 ±0.11, 
Y L /Y = 0.52 ±0.07 ±0.04, (24) 
\P\ 2 = 0.16 ±0.08 ±0.04. 

Evidently, the new data of Yl/Y and R tend to decrease. As a consequence, the comparison 
between theory and experiment is turned the other way around: While several model predic- 
tions for Yl/Y are consistent with CLEO or CDF, almost all models (except for CDDFGN) 

2 The previous CLEO result (27): T L /r = 0.80 ± 0.08 ± 0.05 is based on a subset of the data used in this 
complete analysis. 
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fail to accommodate a small R. It appears that IW(ii) and LF models are most close to the 
data: their results for Tl/T and \P\ 2 are in agreement with experiment, but the predicted R 
is too large by 2 standard deviations compared to the experimental central value. Hence, the 
factorization approach with x Al = Xa 2 = X v can account for the longitudinal polarization 
fraction and parity-odd transverse polarization observed in B — > J/ipK*. In order to explain 
the production ratio and Tl/T simultaneously, it is clear from Eq.(14) that we need f\ 

(25) 



Aa 2 



noting that c\/a2 > 0. Therefore, the data of B — > J/ipK(K*) can be understood provided 
that the nonfactorizable terms in B — > J/ipK* are about the same in S-, P-, and D waves 
and that Xf m B — > J/ipK is slightly larger than that in B — > J/ipK*. The relativistic 
LF quark model should be more reliable and trustworthy than the IW(ii) model since the q 2 
behavior of form factors is directly calculated in the former. We will thus confine ourselves 
to the LF model in ensuing discussion. 

To determine the magnitude of nonfactorizable terms we have to consider the branching 
ratios. The decay rates of B — > J/ipK and B — > J/ipK* are given by 

a 2 G F V cs V* b f J /^F 1 



T(B -> J/i)K) 



47T 



1.102 x 10~ 14 GeV \af F( 



,eff 



" h { m J/i>) 



(26) 



and 



T(B -> J/ipK* 



Pc 



167rmR 



afGpVcsV^jj/^mj/^mB + mj/^)A 



BK* , 
1 1 



rn 



J/i" 



bxf + 2(1 + cV 



(27) 



1.19 x 10~ i4 GeV \af Af K * 



bxf + 2(1 + c 2 y 2 ) 



where p c is the cm. momentum, and uses of fj/^ = 394 MeV extracted from J/ijj — > e + e 



and \V r 



0.038 have been made. Fitting (26) and (27) to the branching ratios 

;i.08± 0.09 ±0.09) x 10~ 3 ; 
1.01 ±0.14) x 10" 3 , 
0.92lg;^±0.08) x 10~ 3 ; 
1.15 ±0.23 ±0.17) x 10^ 3 , 
1.41 ±0.20 ±0.24) x 10" 3 ; 
1.58 ±0.47 ±0.27) x 10~ 3 , 
1.32 ±0.15 ±0.17) x lO" 3 ; 
1.36 ±0.27 ±0.22) x 10~ 3 , 



B(B + - 


-» J/tpK + ) 


B(B° - 


-> J/*PK ) 


B{B + - 


> J/i)K* + ) 


B(B° - 


■> J/ipK* ) 



(2t 



3 For other form-factor models, one can always adjust small nonfactorizable contributions (for example, 
Xa > Xa 2 > Xv) to accommodate the data, as illustrated in [|l0[ [l3|, |9). However, in practice, it is very 
difficult to pin down Xa ± > Xa 2 ano - X v separately if they are not the same. Our results (25) and (31) 
correspond to a minimal modification to the factorization hypothesis. 

4 Since we have taken into account the q 2 dependence of the fine-structure constant, our values of fj/^ 
and f^' given below are slightly larger than the values cited in the literature. 
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where the upper entry refers to the CLEO data || and the lower entry (except for B(B + — > 
J/ijjK + ) which comes from the PDG ||) to CDF data 0, we find averagely 

af(B -> J/ipK) = 0.30 , af(B -> J/ij)K*) = 0.26 , (29) 

for the lifetimes || 

r(B°) = 1.56 x 10~ 12 s, r(5 ± ) = 1.62 x 10~ 12 s, (30) 

and the form factors (21). Using ci(m 6 ) = 1.12 and c 2 {m b ) = -0.27 for Aq^d = 200 MeV, 
we obtain 

X Fl = 0.17, Xai ~ Xa2 ~ Xv = o.U. (31) 

Since |a| ff (S — > J/tpK^)\ is numerically close to |c2(mb)|, one may tempt to argue that 
the large- N c approach works also for B decays. However, this possibilty is ruled out by the 
observed constructive interference in charged B ± — > D^n(p) decays f2"7fl , which implies a 



positive ratio a 2 /ai. Therefore, the sign of a| ff in (29) and X 2 in (31) is fixed to be positive, 
in dramatic contrast to the charm case. The result (31) supports the conjecture (13), but 
the question of why the nonfactorizable term, which lies in the range —0.60 < X2 < — 1/3 
for charm decay fL3| |, becomes positive in B decay remains mysterious and baffling. (For a 



recent attempt of understanding the positive sign of a| ff and X 2 in B decay, see |3( 

With the result (29) or (31) we are ready to discuss B — > ip'K(K*) decays where ip' = 
ip(2S). The analytic expressions for analogus T' L /T', R' and \P'\ 2 are the same as (14), (15) 
and (20) except that the coefficient 1.08 in (14) is replaced by 2.45 and that the unprimed 
quantities a, b, c, x, y, z are replaced by primed ones: 



a = 2.051, b' = 0.733, c' = 0.356, 



We find in the LF model that 

F? K {ml,) = 0.92, AfK'irrfy) = 0.44, Af^mJO = 0.55, V BK \m\ t ) = 0.66. (33) 

Using = 293 MeV determined from ip' — > e + e~ , and assuming a| (B — > ip'K) ~ a| (B — > 
J/ijjK) and af(B -> ip'K*) ~ af(B -> J/ifrK*), we obtain 



.-3 



and 



B(B° -> ^'if ) = 0.50 x KT 3 , £(£+ -> = 0.52 x 10" 

B(B° -> ^*°) = 0.76 x 10~ 3 , 0(5+ -> ^* + ) = 0.79 x 10" 3 , (34) 



r^/r'=0.33, #=1.57, |P'| 2 =0.15. (35) 
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Our results for branching ratios agree with the new CDF measurements [[|: 

B(B + ^ip'K + ) = (0.68 ± 0.10 ± 0.14) x 10" 3 , 

B^^^j'K* ) = (0.90 ±0.21 ±0.20) x 10~ 3 . (36) 

Note that our prediction of T' L /T' is quite different from the predicted range 0.50 < T' L /T' < 



0.67 given in 31 



Finally we trun to the decay mode B® —>■ J/ifxfr. Following [25|, the relevant form factors 
are calculated in the LF model to be 

A^im 2 ^) = 0.35 , A^(m^) = 0.39 , V B °*(m 2 J/f ) = 0.49 . (37) 

A straightforward calculation yields 

T L /T(B S -> J/U) = 0-51 , \P\\B S -> J/M) = 0.13 , (38) 

and 

B(B S -> J/ip(f)) = 1.29 x 10~ 3 (39) 

for ~ 0.25. Experimentally, longitiudinal polarization fraction and the branching ratio 
have been measured by CDF: 

r L /r = o.56 ± o.2i±g:S l 28 l b (b s -»• J/M) = (0.93 ±0.28 ±0.17) x io^ 3 @. (40) 

To summarize, naive factorization or new factorization with universal nonfactorizable 
terms cannot explain the data of Tl/T, R and \P\ 2 simultaneously in all existing form-factor 
models. However, we have shown that a minimal modification to the factorization hypothesis 
can accommodate the recently available data from CDF and CLEO: the nonfactorizable 
term ^ in B — > J/ipK should be slightly larger than x Al ~ Xa 2 ~ X v * n B ~* J/i^K* 
decays. When generalized to B — > ip{2S)K{K*) and B s — > J/ip4>, our predictions agree 
with presently available measurements. Generally speaking, the bulk of exclusive decays 
of heavy mesons can be grossly accounted for by the factorization approach with universal 
nonfactorizable contribution x ^ ~ f f° r charm decays and with x ~ 0.15 for bottom 
decays; that is, the new factorization scheme is very different for charm and bottom decays. 
Just as the charm case, the data of B — > J/if)K(K*) start to reveal departures from the 
factorization postulation: the nonfactorizable term x m B decay is also process dependent 
and not universal. 
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